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^ Abstract 
ON 

Predictions of elliptic flow (t^) and nuclear modification factor (Raa) are provided as a function of 

^ centrality in U + U collisions at y/s NN = 200 GeV. Since the 238 U nucleus is naturally deformed, one 

' | could adjust the properties of the fireball, density and duration of the hot and dense system, for example, 

O in high energy nuclear collisions by carefully selecting the colliding geometry. Within our Monte Carlo 

£H Glauber based approach, the v 2 with respect to the reaction plane vf p in U + U collisions is consistent 

1 with that in Au + Au collisions, while the v 2 with respect to the participant plane v 2 p increases ~30- 

{N| 60% at top 10% centrality which is attributed to the larger participant eccentricity at most central 

+ U collisions. The suppression of Raa increases and reaches ~0.1 at most central U + U collisions 

O that is by a factor of 2 more suppression compared to the central Au + Au collisions due to large size 

^ and deformation of Uranium nucleus. 
O 
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^ 1. Introduction 

X ■ n 

Most striking findings at RHIC are the large elliptic flow v 2 [l| and the strong suppression of nuclear 
. .^modification factor Raa 12] • The v 2 is defined by the second harmonic Fourier coefficient of azimuthal 
particle distribution with respect to the reaction plane, and the Raa is defined by the ratio of invariant 
yield in A + A collisions to that in p + p collisions scaled by number of collisions. Recent systematic 
measurements of v 2 ^ as well as developments of viscous hydro dynamical models [H, H, @, 0] provide 
a conservative upper limit of the viscosity rj to the entropy s ratio r]/s < 0.5. This corresponds to 
the 6 times larger value of an absolute lower bound 77/s = l/An predicted by strongly coupled gauge 
filed theories based on the AdS/CFT correspondence @, 0]. It has been observed that the ratio v 2 /s 
in different systems from AGS to RHIC scale like 1/ SdN^/dy [l(| as it was predicted by a low density 



limit of t>2 llll; ll2|]; where e is the initial geometrical anisotropy (eccentricity), S is the transverse area 
and dN ch /dy is the charged particle rapidity density. The saturation of v 2 /e would indicate that the 
system is approaching the hydrodynamical limit and the collectivity no longer increases when the system 
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size becomes larger. The measurements of transverse momentum spectra of charged hadrons showed 
that the yield at most central Au + Au collisions at ^s NN = 200 GeV is suppressed by a factor of ~5 



compared to the p + p reference scaled by number of binary collisions iV co ii [13l . [I4j and the similar 
level of suppression persists for neutral pions up to px = 20 GeV/c [l5j]. The integrated Raa above 
Pt > 5 GeV/c and > 10 GeV/c decreased monotonically as a function N pax t and there were no sign of 



saturation 15 



Assuming the underlying dynamics remains the same, we ask what would happen to v 2 an d Raa 
for a larger colliding system 238 U + 238 U collisions ? Comparing to the 197 Au nucleus, the 238 U has a 
much larger mass and, more importantly, it is largely deformed. The planned U + U collisions at RHIC 
will be important for us to understand how those observables behave at higher particle density. Monte 
Carlo Glauber simulations showed that the transverse number density 1/ SdN^/dy increases ~ 35% at 
most central events in ideal tip-tip collisions (head-on collisions along the longest axes) (l6|^The U + U 



collisions will become possible when the Beam Ion Source becoming operational in 2012 |171 . 

In this letter, we will report a geometrical approach based on the Monte Carlo Glauber model to 
predict the elliptic flow v 2 as well as the nuclear modification factor Raa in U + U collisions at top RHIC 
energy. In the Section [2J we will discuss our parameterization of Glauber model and define geometrical 
quantities which are used in this study. In the Section El the results of v 2 and Raa in U + U collisions 
will be presented and compared to the data in Au + Au collisions at y/s NN = 200 GeV. 



2. Glauber Model 



The nucleon density distribution is parameterized by a deformed Woods-Saxon profile [18 
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where po is the normal nuclear density, R and a denote the radius of nucleus and the surface diffuseness 
parameter, respectively. We have used R = 6.38 fm and a = 0.535 fm for Au nucleus, and R = 6.81 fm 
and a = 0.55 fm for 238 U nucleus. The Y™(9) denotes the spherical harmonics and 9 is the polar angle 



with the symmetry axis of the nucleus. Deformation parameters are (3 2 = 0.28 [19[ and (5^ = 0.093 [20 
for Uranium. The presence of modifies the shape of Uranium compared to that only with /3 2 , 

The radius increases ~6% (3%) at 9 = 



which was implemented in several different models [16|, [19 
{9 = 7r/2), while it decreases ~3% around 9 = it /A. We have assumed that Au nucleus is spherical 
(f3 2 = 04 = 0), thus Eq. ([T]) reduces the spherical Woods-Saxon profile. Recent calculation [2l| shows 
that the ground-state deformation of 197 Au affects the eccentricity of initial geometry overlap only 
at most central collisions from both optical and Monte Carlo Glauber simulations. The positions of 
nucleons are sampled by 4nr 2 sin (9)p(r) d9d<j), where the absolute normalization of p(r) is irrelevant. 

Both projectile and target U nuclei are randomly rotated along the polar and azimuthal directions 
event-by-event with the probability distribution sin and uniform distribution for G and $, respectively. 
The sin G weight needs to be implemented to simulate unpolarized nucleus-nucleus collisions. The results 
are averaged over all possible orientations unless otherwise specified. 
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A binary nucleon-nucleon collision take places if 



d < 



O'NN 



71 



(3) 



where d is the distance between nucleons in the transverse direction orthogonal to the beam axis, and 
<jnn = 42 mb is the inelastic nucleon-nucleon cross section at y/s = 200 GeV. For each event, the total 
number of binary collisions iV co ii is calculated by the sum of individual number of collision and the total 
number of participant nucleons N paxt is the number of nucleons that interacts at least once. 



Charged particle pseudorapidity density is obtained by a two component model [22 
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where n pp = 2.29 and x = 0.145 are fixed to reproduce the PHOBOS results [23j. Event-by-event 
multiplicity fluctuations have been taken into account by convoluting Negative Binomial Distribution 
for a given iV part and 7V coU 
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where \x 



n 



pp 



is the mean of the distribution and 1/k = 0.5 corresponds to deviation from a Poisson 



distribution. In this study, we have generated 1 million events for U + U collisions by randomly selecting 
an impact parameter b according to the do/db = 27rb. 

Figure [1] shows the comparison of dN^/drj distributions in Au + Au and U + U collisions at y/s NN 
= 200 GeV from our Monte Carlo Glauber model. The maximum dN^/dn in U + U collisions increases 
~15% compared to that in Au + Au collisions. We have defined the event centrality bins by the fraction 
of events in dN^/dn. The centrality bins are summarized in Table [TJ 



Tabic 1: Summary of centrality bins based on the dN^/dr] and (N pai t), (N co \\), (Sr.p), (Spp), (erp), £pp{2} = vv^pp) 
and (L) for each centrality bin in U + U collisions at y/s NN = 200 GeV. Errors denote systematic uncertainties, see texts 
for more details of systematic error evaluations. 



centrality dN ch /drj (N paxt ) 


(iV coll ) 


(Srp) 


(5 PP ) 




epp{2} 


(L) (fin) 


0-5% 


> 740 


418 ± 6 


1341 ± 105 


30.9 ± 1.7 


29.7 ± 1.7 


0.021 ± 0.007 


0.156 ± 0.004 


4.4 ± 0.1 


5-10% 


> 609 


358 ± 14 


1058 ± 52 


27.1 ± 1.9 


26.9 ± 1.9 


0.08 ± 0.02 


0.18 ± 0.01 


4.2 ± 0.2 


10-20% 


> 410 


281 ± 13 


751 ± 49 


22.9 ± 1.8 


22.7 ± 1.8 


0.15 ± 0.02 


0.24 ± 0.02 


3.9 ± 0.1 


20-30% 


> 269 


199 ± 14 


462 ± 45 


18.4 ± 1.6 


18.2 ± 1.6 


0.23 ± 0.03 


0.31 ± 0.03 


3.5 ± 0.1 


30-40% 


> 170 


137 ± 14 


272 ± 39 


14.8 ± 1.5 


14.5 ± 1.6 


0.29 ± 0.03 


0.38 ± 0.03 


3.2 ± 0.2 


40-50% 


> 101 


89 ± 13 


149 ± 31 


11.8 ± 1.5 


11.4 ± 1.5 


0.34 ± 0.04 


0.45 ± 0.04 


2.9 ± 0.2 


50-60% 


> 56 


55 ± 11 


75 ± 22 


9.3 ± 1.5 


8.8 ± 1.5 


0.38 ± 0.04 


0.51 ± 0.05 


2.6 ± 0.2 


60-70% 


> 29 


31 ± 9 


35 ± 13 


7.1 ± 1.5 


6.5 ± 1.6 


0.39 ± 0.05 


0.59 ± 0.07 


2.3 ± 0.2 


70-80% 


> 13 


16 ± 6 


15 ± 8 


4.8 ± 1.7 


4.0 ± 1.8 


0.38 ± 0.06 


0.68 ± 0.09 


1.9 ± 0.3 



3 




dNJdn 



Figure 1: (Color online) dN^/dr] distribution in Au + Au (solid line) and U + U collisions (filled circles) at y/s NN = 200 
GeV by averaging over all orientations. 
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Since the positions of nucleons fluctuate event-by-event, the principal axes of the participant nucleons 
in the transverse plane are tilted and rotated with respect to the original coordinate system. We define 
the participant plane (PP) which is the relevant plane to take into account the event-by-event position 
fluctuations of participant nucleons. The transverse area and eccentricity with respect to the reaction 
plane (RP) and participant plane are defined as 



Srp = nJ<T*<T* (6) 



Spp = 7iy a%a 2 - a 2 y , (7) 
£rp 




where a 2 = {x 2 } — {x} 2 , a 2 = {y 2 } — {y} 2 and a xy = {xy} — {x}{y}. The curly brackets {...} denote 
the average over all participants for a given event. We have also calculated the averaged transverse path 
length L from the RMS width 

'^ + ^, (10) 

which could be a relevant geometrical quantity for the Raa- The path length is very close to pL defined 
in 24j, while Eq. ({TO]) takes into account the event-by-event center of mass shift of the nuclei within 
the transverse plane. Average quantities, (N paxt ), (N co \\), (S), (e) and (L) have been calculated for each 
centrality bin where (...) describe the average over all events. 

Systematic uncertainties on the average quantities have been estimated (i) by varying input param- 
eters R, well as the total cross section within ±5% and (ii) by using different density profiles 
for nucleons in the Monte Carlo Glauber simulations. The dominant source of systematic uncertainty 
is the total cross section. Total systematic uncertainty has been evaluated by the quadratic sum of 
individual systematic uncertainty. Table Q] summarizes the centrality bins, average quantities and their 
systematic uncertainties obtained in the Monte Carlo Glauber simulation. 

Figure [2] compares the iV co n, transverse number density 1/ (Snp)dN c i l / 'drj, reaction plane eccentricity 
(£rp), second order cumulant of participant eccentricity epp {2} = v\^pp)> an d path length (L) as a 
function of N part together with their systematic uncertainties in Au + Au and U + U collisions at y/s NN 
= 200 GeV. We found that the all geometrical quantities essentially scale with iVp art . The iV co n, (S), (L) 
increase and (srp) decreases at most central U + U collisions compared to those in Au + Au collisions 
because of the larger size of Uranium. One can see that the £pp{2} in U + U collisions starts deviating 
the iVp ar t scaling around iVp art = 200, and increases ~60% at top 5% central U + U collisions. The 
higher values of £pp{2} in U + U collisions for large N pait is purely from the ground-state deformation 
of Uranium. We have confirmed that the £pp{2} becomes the same if we assume the Uranium is 
spherical. The relevance of the epp {2} will be discussed in the next section. 
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Figure 2: (Color online) Comparison of (a) N co n, (b) l/(S-&p)dN c i 1 /dr], (c) (erp), (d) epp{2} and (c) (L) as a function of 
-Wpart in Au + Au (open circles) and U + U collisions (filled circles) at y/s NN — 200 GeV. Open boxes show the systematic 
uncertainties in Monte Carlo Glauber simulations. 
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3. Results and Discussions 



3.1. Elliptic flow t>2 

It has been found that the elliptic flow v 2 divided by initial anisotropy e in coordinate space scaled 
like 1/ SdN c h/dy among different energies and collision systems from AGS to RHIC [loj. A simple 



formula that has been proposed in (25| describes very well the variation of V2 with 1/SdN/dy 

v 2 _ h 

7 ~ 1 + B (l/SdN/dy)- v 



where dN/dy is the rapidity density of total particles, h is the V2/S in the ideal hydrodynamical limit 
when 1/S dN/dy — > 00, and B contains informations about the equation of state and the partonic cross 
section [25[. The Eq. (TTTT) reduces v 2 /e ~ (h/B)l/S dN/dy when 1/S" dN/dy — > for leading order 
in 1/S dN/dy, thus the above equation satisfies both low density and ideal hydrodynamical limit of 
i>2- The integrated V2 for unidentified charged hadrons from the PHOBOS collaboration can be well 



described by the Eq. ( TTTT) 26]. Assuming no change in the collision dynamics, we will study the i^/e 
distributions versus the collision centrality in U + U collisions. 

Figure [3] shows the v 2 {4:} / (e-rp) and v 2 {EP} / £pp{2} as a function of 1/ (S)dN/dy in Au + Au 
collisions at x /s NN = 200 GeV, where v 2 {4:} and t^{EP} denote the V2 from four particle cumulant 



method [27|, |28| and that from standard event plane method [29| , respectively. The dN/ dy is obtained 
by multiplying 3/2 to the measured dN^/dy at STAR to take into account the neutral particles. 
The simultaneous fit has been performed for (a) four particle cumulant t>2{4}, six particle cumulant and 
g-distribution method and for (b) standard event plane ^{EP}, scalar product method, rj subevent, 



random subevent, and two particle cumulant. The results of v-i are taken from [3JJ. The two different 
groups of V2 are categorized based on the multi-particle methods for (a) and two particle methods for 
(b). As long as the distribution of eccentricity is 2D gaussian in the transverse plane, the effect of fluc- 
tuation on the v 2 {4:} is negligible and thus the £rp can be used to scale the ^{4} j32[. This assumption 
holds except for the peripheral 60-80% centrality, where the distribution of eccentricity becomes non- 
Gaussian. The V2 from two particle methods are expressed as (v% where t>2 is the true i> 2 value and 

a varies from 1 to 2 depending on the event plane resolution (ijl, 34]. In this study, £ par t{2} = yj (£ part ) 

was used by assuming a = 2. We confirmed that the resulting v 2 values unchanged by using e part (i.e. 
a = 1). Because the v 2 values were extrapolated from fa/e) multiplied by e, most of the difference 
between e pa rt{2} and 5 part is canceled out and thus the resulting w 2 is the same regardless of the choice 
of eccentricity. 

Figure H] shows the extracted v-i in U + U collisions compared to those in Au + Au collisions as a 
function of centrality at y/s NN = 200 GeV. The vf p (v p p ) denotes the V2 measured with respect to the 
reaction (participant) plane. The vf p and v PF have been calculated by multiplying the (£rp) and 5pp{2} 
to the fitting results of v^/s shown in Fig.[3]for each centrality bin. Since we have calculated the dN^/drj 
in the Monte Carlo Glauber simulation, it is necessary to convert the dN^/drj to dN/dy for calculating 
the i>2 for each centrality bin. We assume that dN^/dy ~ lAbdN^/drj to extrapolate the v^/e for each 



centrality 35|. The vf p in U + U collisions is consistent with that in Au + Au collisions for centrality 
0-80%. The v P p is also consistent with each other in U + U and Au + Au collisions for centrality 20- 
80%, whereas the v 2 in U + U collisions at top 0-10% centrality is 30-60% larger than that in Au + Au 
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Figure 3: (Color online) (a) Four-particle cumulant V2{4}/(erp) as a function of l/(Spp)dN/dy for unidentified charged 
hadrons in Au + Au collisions at y/s^ = 200 GeV. (b) The same plot as (a) for the standard event plane u 2 {EP}/epp{2} 
as a function of 1/ (Spp)dN/dy. Both the values of t>2{4} and 1>2{EP} are taken from [3lj]. Only statistical errors on the 
V2 are shown and are smaller than symbols. Open boxes are systematic errors from the Monte Carlo Glauber simulation. 
Solid lines are fitting results by Eq. (fTTj) . See more details about fitting in the texts. 
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Figure 4: (Color online) (a) The i;^ p as a function of centrality. Solid circles is the t>2{4}, open boxes and shaded band 
show the extracted v-i from the fitting of V2/(£rp) and (erp) in Monte Carlo Glauber simulation in Au + Au and U + U 
collisions, respectively, (b) The same comparison of i> PP with £pp{2} as (a), where solid stars are the ^{EP}. The errors 
on the V2 include systematic errors from Monte Carlo Glauber simulation and errors from fitting of vijs. 



9 



collisions. The larger v 2 is attributed to the larger participant eccentricity due to the ground-state 
deformation in top 0-10% centrality in U + U collisions compared to that in Au + Au collisions. The 
extracted t>2 in Au + Au collisions are slightly smaller than the data at peripheral collisions. Since the 
dN^/drj has been tuned to reproduce the PHOBOS results and is smaller than the STAR dN^/dy at 
peripheral 60-80%, the resulting V2/e (and hence the v%) become smaller than the STAR vi- 

3. 2. Nuclear modification factor Raa 

The integrated Raa over a certain p T range in Au + Au collisions at y/s NN = 200 GeV has been 
described by Raa = (1 — SoN pSuTt a ) n ~ 2 , where n — 8.1 is the power-law exponent of pt distribution, and 
S = (9.0±6.1) x 10~ 3 and a = 0.57±0.14 for A part > 20 andp T > 5 GeV/c 0. We have assumed that 
the path length (L) determines the Raa in both Au + Au and U + U collisions. The Raa in U + U 
collisions has been extrapolated by fitting the Raa{L) in Au + Au collisions with an ansatz from above 
equation 

R AA (L) = (1 - S' {L) b r~ 2 , (12) 
where n = 8.1, S' and b are free parameters that have been evaluated by fitting the data. 



Figure [5] shows integrated Raa for > 5 GeV/c from [15j as a function of (L) in Au + Au collisions 
a ^ \/ s nn = 200 GeV. We have assumed that the definition of our centrality bins is the same as that of 
the PHENIX in order to plot the Raa as a function of (L) for each centrality. Result of the fit with 
Eq. Q12p is shown by the solid line in Fig. [5] and holds quite well over entire range of (L) since the (L) 
scales like Apart 1 / 3 . 

Figure [6] shows integrated Raa as a function of A part extrapolated for U + U collisions at y/s NN 
= 200 GeV. The Raa in U + U collisions has been evaluated for a given (L) in each centrality. The 
calculated Raa i n Au + Au collisions (open boxes) is consistent with the data within the systematic 
error as it should. We found that the Raa reaches ~0.1 at most central U + U collisions, which is by 
a factor of 2 more suppression compared to the central Au + Au collisions due to the larger size and 



the deformation of Uranium. Heinz and Kuhlman pointed out in [19| that the radiative energy loss AE 
of a fast parton moving through the medium is almost independent of the orientations of nuclei for the 
out-of-plane direction in the full overlap U + U collisions. Whereas the AE in the in-plane direction 
decreases by about 35% towards the ideal body-body collisions (head-on collisions along the shortest 
axes). For the body-body collisions, they found that the difference of AE between out-of-plane and 
in-plane directions is more than twice in U + U collisions that achieved in Au + Au collisions. They 
also found that the total energy loss is larger by up to a factor of 2. More differential study, such as 
selecting the orientations of Uranium and directions with respect to the reaction plane, will be needed 
to see whether the Raa would have such dependences or not. Since the (L) in U + U collisions is 
slightly larger (~3% in central, and ~5% in peripheral collisions) than that in Au + Au collisions, the 
Raa in U + U collisions would be even more suppressed for a given A part . Due to the large errors on 
the extrapolated Raa, we have not observed any difference of the Raa between Au + Au and U + U 
collisions for a given A part . 

4. Summary 

In summary, we have predicted the v 2 and Raa in U + U collisions at y/s NN = 200 GeV by a simple 
geometrical approach with the Monte Carlo Glauber simulation. We found that the v^ p is consistent 
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Figure 5: (Color online) Integrated Raa above > 5 GeV/c [l5j as a function of (L) in Au + Au collisions at y/s NN 
= 200 GeV. Most central bin denote 0-10% rather than 0-5%. Statistical and px-uncorrelated errors are smaller than 
symbols. Open boxes are the errors on (L) in x axis and px-correlated systematic errors on the Raa m y axis. The 
dashed lines and the single box on right at unity show the errors on A^ co n and normalization of the p + p reference, 
respectively [Hj]. Solid line is the fitting result by Eq. (fT2|) . 
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circles) [HI are plotted for comparison. 
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with that in Au + Au collisions over all centrality range, whereas the v% p increase by 30-60% at most 
central 0-10% collisions due to the larger e PP {2} in U + U collisions. The Raa at top 5% central U + U 
collisions further suppressed and reaches ~0.1, which is by a factor of 2 more suppression compared 
to the most central Au + Au collisions. It is clear that the larger mass and deformation form the U 
nucleus will allow us to study the matter at higher density. By selecting the relative orientation of the 
colliding Uranium nuclei, the discussed effects may be further enhanced. We will report the method in 
a separate paper. 
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